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Abstract 



An effective Hamiltonian consisting of bare A «-> ttN, jN vertex interactions 

p ' 

and energy-independent meson-exchange ttN *-* ttN, jN transition operators 

is derived by applying a unitary transformation to a model Lagrangian with 

Qs, ' N,A,tt, p, to, and 7 fields. With appropraite phenomenological form factors 
O ' 

© ; and coupling constants for P and A, the model can give a good description 

j of ttN scattering phase shifts up to the A excitation energy region. It is 

On" 

shown that the best reproduction of the recent LEGS data of the photon- 
asymmetry ratios in 7p — ► ir°p reactions provides rather restricted constraints 
on the coupling strengths Ge of the electric E2 and Gm of the magnetic 
Ml transitions of the bare A <-> 7^ vertex and the less well-determined 

X" 

I— ( ' coupling constant g^NN of u meson. Within the ranges that Gm = 1.9 ± 



0.05, G E = 0.0 ± 0.025, and 7 < g uNN < 10.5, the predicted differential 
cross sections and photon-asymmetry ratios are in an overall good agreement 
with the data of 7^ — ► 7r°p, jp — > 7r + n, and jn — > ir~p reactions from 180 
MeV to the A excitation region. The predicted M 1 + and E<±+ multipole 
amplitudes are also in good agreement with the empirical values determined 
by the amplitude analyses. The constructed effective Hamiltonian is free of 
the nucleon renormlization problem and hence is suitable for nuclear many- 
body calculations. We have also shown that the assumptions made in the 
iT-matrix method, commonly used in extracting empirically the 7^ — > A 
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transition amplitudes from the data, are consistent with our meson-exchange 
dynamical model. It is found that the helicity amplitudes calculated from 
our bare 7iV — > A vertex are in good agreement with the predictions of the 
constituent quark model. The differences between these bare amplitudes and 
the dressed amplitudes, which are closer to the empirical values listed by 
the Particle Data Group, are shown to be due to the non-resonant meson- 
exchange mechanisms. Within the range 7 < g^NN < 10.5 of the oj meson 
coupling favored by the data of the photon-asymmetry ratios in jp — > n°p 
reactions, our values of the E2/M1 ratio for the 7./V — > A transition are (0.0 
± 1.3)% for the bare vertex and (-1.8 ± 0.9)% for the dressed vertex. 



PACS numbers: 13.75.Gx, 24.10.1, 25.20.Lj 
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I. INTRODUCTION 



The main objective of investigating photo- and electro-production of mesons on the 
nucleon is to study the structure of the nucleon excited states (N*). This has been pursued 
actively |l| during the period around 1970. With the developments at several electron 
facilities since 1980, more extensive investigations of the A-excitation have been carried 
out both experimentally and theoretically 0. Apart from the need of precise and extensive 
measurements which will soon be possible at CEBAF and Mainz, an accurate understanding 
of the N* structure can be obtained only when an appropriate reaction theory is developed 
to separate the reaction mechanisms from the hadron structure in the 'jN — > ttN, ttiiN 
reactions. The importance of this theoretical effort can be understood by recalling many 
years experiences in the development of nuclear physics. For example, the information 
about the deformation of 12 C can be extracted from 12 C(p,p') 12 C*(2 + , 2.44 MeV) inelastic 
scattering only when a reliable reaction theory j3], such as the Distorted- wave Impulse 
Approximation or the coupled- channel method, is used to calculate the initial and final 
proton- 12 C interactions. Accordingly, one expects that the N* structure can be determined 
only when the interactions in its decay channels jN, ttN, and 7T7riV can be calculated 
from a reliable reaction theory. It is the objective of this work to address this problem 
from the point of view of meson-exchange models. In contrast to approaches based on the 
dispersion relations |J or the .fT-matrix method [§-0] , our approach is aimed at not only an 
investigation of the N* structure but also on the application of the constructed model to a 
consistent calculation of N* in nuclear many-body systems. 

The meson-exchange models have been very successful in describing nucleon-nucleon 
interactions J8|, electroweak interaction currents [P,flO[|, meson-meson scattering and 



meson-nucleon scattering |12|- |14T| . It is therefore reasonable to expect that the same success 
can also be achieved in the investigation of pion photo- and electro-production. This possi- 
bility has, however, not been fully explored. The dynamical models of pion photoproduction 
developed in Refs. [|l5|-|T7|j did contain the well-established meson-exchange mechanisms of 
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pion photoproduction, but phenomenological separable potentials were used to describe the 
ttN multiple scattering. The improvement made in Ref. ]Tj| suffered from the theoretical 
inconsistency in defining the meson-exchange irN interaction and ttN — > 7iV transition. 



The model developed in Refs. flT9"l , |20f also does not treat meson-exchange completely since 
a zero-range contact term is introduced to replace the particle-exchange terms of their ttN 
potential. In all of these models, the incomplete treatment of the meson-exchange interac- 
tions leads to some uncertainties in interpreting the parameters characterizing the 'jN — > A 
vertex which is the main interest in testing hadron models. The formulation developed in 
Ref. |2]J can, in principle, be used to examine the meson-exchange mechanisms in pion pho- 
toproduction, but has not been pursued numerically. In this work, we will try to improve the 
situation by appling the unitary transformation method developed in Ref. to derive from 
a model Lagrangian an effective Hamiltonian for a consistent meson-exchange description 
of both the ttN scattering and pion photoproduction. Furthermore, the constructed model 



can be directly used to improve and extend the irNN Hamiltonian developed in Ref. |23| to 
also describe the electromagnetic A excitation in intermediate energy nuclear reactions. 

The starting point of constructing a meson-exchange model is a model Lagrangian of 
relativistic quantum field theory. The form of the Lagrangian is constrained by the observed 
symmetries of fundemental interactions, such as Lorentz invariance, isospin conservation, 



chiral symmetry, and gauge invariance. The most common approach ]24j is to find an ap 



propriate three-dimensional reduction of the Ladder Bethe-Salpeter equation of the consid- 
ered model Lagrangian. The meson-exchange potentials are then identified with the driving 
terms of the resulting three-dimensional scattering equation. The most recent examples are 



the irN models developed in Refs. [12-14, 19|. The extension of this approach to investigate 



pion photoproduction has also been made in Ref. |20[ . 

Alternatively, one can construct a meson-exchange model by deriving an effective Hamil- 
tonian from the considered model Lagrangian. Historically, two approaches have been de- 
veloped. The first one is to use the Tamm-Dancoff approximation [|25|| . This method leads 
to an effective Hamiltonian which is energy-dependent and contains unlinked terms, and 



hence can not be easily used in nuclear many-body calculations. A more tractable approach 
is to apply the method of unitary transformation which was developed by Fukuda, Sawada 
and Taketani P6| , and independly by Okubo [P? |. This approach, called the FST-Okubo 



method, has been very useful in investigating nuclear electromagnetic currents p8|-|30|j and 
relativistic descriptions of nuclear interactions pTLpi]. The advantage of this approach is 
that the resulting effective Hamiltonian is energy independent and can readily be used in nu- 
clear many-body calculation. Motivated by the investigation of the nNN dynamics p3| , |35 



this method has been extended in Ref. |22| to derive an effective theory involving pion pro- 



duction channels. In this work, we adopt this method to develop a dynamical model for ttN 
scattering and 'yN — > ttN reactions. 

It is necessary to explain here how our approach is related to the approach based on 
the Chiral Perturbation Theory (CHPT) |p6| . Since chiral symmetry is a well-established 
dynamical symmetry of strong interactions, it should be used to constrain our starting 
Lagrangian. This leads us to assume that our starting Lagrangian is an effective Lagrangian 
for generating the tree-diagrams in CHPT. The parameters are then completely determined 
by the well-established chiral dynamics such as PCAC and current algebra. Therefore, 
our model and CHPT are identical in leading orders. The differences come from how the 
unitarity is implemented to account for the irN multiple scattering. In the spirit of CHPT, 
the "low" momentum pions are considered as weakly interacting Goldstone bosons and hence 
their interactions with the nucleon can be treated as perturbations ||37|| . This amounts to 
restoring the unitarity perturbatively by calculating loop corrections order by order. It 
is then necessary to include more terms in the effective Lagrangian. A phenomenological 
procedure is then unavoidable to determine the accompanied low-energy constants. 

In the meson-exchange model, one hopes to describe the ttN multiple scattering in the 
entire kinematic region including the highly non-perturbative A excitation region. The es- 
sential assumption is that the 7rN multiple scattering is governed by a few-body Schrodinger 
equation with the driving terms calculated from the starting Lagrangian in a perturbation 
expansion in the coupling constants. This can be realized in practice only when the driving 



terms are regularized by appropriate phenomenological form factors. Qualitatively speaking, 
the meson-exchange model is an alternative to CHPT in the kinematic region where pertur- 
bative calculations become very difficult or impossible. Both approaches involve phenomeno- 
logical parameters. The success of each approach depends on whether these parameters can 
be interpreted theoretically. 

In this work we will focus on the A excitation and will limit our investigation to the energy 
region where 2n production is negligibly small. By applying the unitary transformation of 



Ref. |22[ to a model Lagrangian for N, A, ir, p, uj, and 7 fields, we have obtained an effective 
Hamiltonian consisting of bare A <-» ttN, 'yN vertex interactions and energy-independent 
ttN <-> ttN, syN transition operators. The ttN scattering phase shifts |38|-fl0|j are used to 



determine the hadronic part of the constructed effective Hamiltonian which has only seven 
parameters for defining the vertices of the meson-exchange 7rN potential and the A <-> ttN 
transition. The strong vertex functions in the 'jN — > ttN transition operator are then also 
fixed. This is a significant improvement over the previous dynamical models |l5|-|lTf in which 
the employed separable potentials have no dynamical relation with the pion photoproduction 
operator. A consistent description of the ttN scattering and 'yN — > ttN transition is crucial 
for separating the reaction mechanisms due to meson-exchange non-resonant interactions 
from the total 'yN — > A transition. 

Once the hadronic part of the effective Hamiltonian is determined, the resulting pion 
photoproduction amplitude has only three adjustable parameters: Gm of magnetic Ml and 
Ge of electric E2 transitions of the bare 7iV — > A vertex, and the less well-determeined 
uNN coupling constant. We will determine these three parameters by considering the most 
recent LEGS data [[II]] of the photon-asymmetry ratios in 7^ — > 7r°p reactions. The resulting 



parameters are then tested against very extensive data in Refs. |42|-f44||. 

It is customary to test hadron models by comparing the theoretical predictions of N* — > 
'jN transition amplitudes with the empirical values listed by the Particle Data Group [PDG] 
|£5fl . Since the first systematic calculation based on the constituent quark model was 



performed, it has been observed that the predicted A — > ^yN transition amplitudies [^-|50 



are significantly smaller than the empirical values listed by PDG | 45| . While the problem 
may be due to the limitations of the constituent quark model, it is necessary to recognize that 
the empirical values of PDG are obtained by using the .ff -matrix method or dispersion 
relation Both approaches contain assumptions about the non- resonant contributions to 
the 7iV — > A transition and must be justified from a dynamical point of view. Within our 
dynamical model, we will address this point concerning the .fT-matrix method. This leads 
us to identify our bare ^yN — > A vertex with the constituent quark model. The dispersion 
relation approach |l],|5lH is defined in a very different theoretical framework and therefore is 
beyond the scope of this investigation. 

In section II, we will use a simple model Lagragian to explain how an effective Hamil- 



tonian can be constructed by using the unitary transformation method of Ref. p2| . The 
method is then applied to realistic Lagrangians to derive in sections III and IV an effective 
Hamiltonian for ttN scattering and pion photoproduction. The equations for calculating the 
tcN scattering and 'yN — > ttN amplitudes are also presented there. The relationship with 
the .fT-matrix method are then established. Results and discussions are given in section V. 
The conclusions and discussions of future studies are given in section VI. 



II. METHOD OF UNITARY TRANSFORMATION 

To explain the unitary transformation method of Ref. |22] (will be referred to as the SKO 
method), it is sufficient to consider a simple system consisting of only neutral pions and 
fictituous cr mesons. The objective is to derive an effective Hamiltonian from the following 
Lagrangian density 

L{x) =L {x)+Li{x), (2.1) 
where Lq(x) is the usual noninteracting Lagrangian, and the interaction term is taken to be 

Lj(x) = -g a7 ^<t>l{x)<t>a{x)- (2-2) 
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The Hamiltonian can be derived from Eq. (2.1) by using the standard method of canonical 
quantization. In the second-quantization form, we obtain (in the convention of Bjorkin and 



Drell f52|]) 

H = H + H It (2.3) 



with 



H = J dk(E n (k)a\k)a(k) + E a (k)tf(k)b(k)), (2.4) 

f dk\dk 2 dk 
Hi = g an7T I 



(2n)^8E 7T (k 1 )E n (k 2 )E a (k) 
x ([2o t (j?i)o(jfe 2 )6(^)^i -k 2 -k) + a^jzi^fybfysfa + k 2 -k) 
+ a{k 1 )a{k 2 )b{k)8{k 1 + k 2 + k)} + [h.c.]), (2.5) 

where a){k) and b\k) are, respectively, the creation operators for tc and a particles, 



E a (k) = ym^ + k 2 is the free energy for the particle a, and [h.c] means taking the hermi- 
tian conjugate of the first term in the equation. We further assume that the mass of the a 
meson is heavier than two-pion mass; i.e. m a > 1m v . 

Because of the intrinsic many-body problem associated with the starting quantum field 
theory, it is not possible to solve exactly the equation of motion for meson-meson scattering 
defined by the above Hamiltonian. A simplification is obtained by assuming that in the 
low and intermediate energy regions, only "few-body" states are active and must be treated 
explicitly. The effects due to "many-body" states are absorbed in effective interaction op- 
erators which can be calculated in a perturbation expansion in coupling constants. This 
few-body approach to field theory was pioneered by Amado |p3[ . In the SKO approach, this 
is achieved by first decomposing the interaction Hamiltonian Hj Eq. (2.5) into two parts 

Hi = Hf + Hf (2.6) 

Hf = I dk ^ dk (a^k^ihWkWh + k 2 -k) + [h.c.]) (2.7) 

J ^/(2n) 3 8E n (k 1 )E 7T (k 2 )E a (k) 

q r dk\dk 2 dk 

J .l(27r)^8E n (k 1 )E 7T (k 2 )E a (k) 



( [2a t (fc 1 )a(fc 2 )&(fc)5(^i -fa-k) + a(fci)a(fc 2 )&(fc)£(fci + fa + k)} + [h.c.]) (2.8) 

The elementary processes induced by Hf are illustrated in the upper half of Fig. I. For 
m a > 2m v , the a — > 7r7r decay and 7T7t — > cr annihilation are 'real processes' and can 
take place in free space. On the other hand, the processes n <-> na and vacumm <-> 7T7ro" 
induced by if^ are 'virtual processes' (lower part of Fig. 1). They can not occur in free 
space because of the energy-momentum conservation. The essence of the SKO method is 
to systematically eliminate the virtual processes from the considered Hamiltonian by using 
unitary transformations. As a result the effects of 'virtual proceesses' are included as effective 
operators in the transformed Hamiltonian. 
The transformed Hamiltonian is defined as 

H' = UHU + (2.9) 
= U(H + Hf + Hf)U + , 

where U = exp(-iS) is a unitary operator defined by a hermitian operator S. By expanding 
U = 1 — iS + ... , the transformed Hamiltonian can be written as 



H' = H + Hf + Hf + [H , iS ] (2.10) 

_l_ ... 



+ [H I ,iS] + ^[[H Q ,iS],iS 



To eliminate from Eq. (2.10) the virtual processes which are of first-order in the coupling 
constant g an7T , the SKO method imposes the condition that 

H? + [H o ,iS] = 0. (2.11) 

Since H is a diagonal operator in Fock-space , Eq. (2.11) clearly implies that iS must have 
the same operator structure of Hf. 

To simplify the presentation, we write H Q as 



Hf = E / F n On dfadfa dk , (2.12) 
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where O n denotes the part containing creation and annihilation operators, and F n is the rest 
of nth term in Eq. (2.8). In the form of Eq. (2.12), the solution of Eq. (2.11) can be written 

as 



iS 



J SnOndhdhdk. (2.13) 



Our task is to find S n by solving Eq. (2.11). Considering two eigenstates \i > and \f > of 
the free Hamiltonian H such that < f\O n \i >= 1, Eqs. (2.11), (2.12) and (2.13) then lead 
to 

iS n = —^—F n . (2.14) 
h f - hi 

Note that Ei and Ef are the eigenvalues of free Hamiltonian H , and hence the solution S n 
is independent of the collision energy E of the total Hamiltonian H. This is an important 
feature distingushing our approach from the Tamm-Dancoff approximation. By using the 
above relation, it is easy to verify that the solution of the operator equation (2.11) is 

dk\dk 2 dk 



iS Qairir j 



^/(2Tr) 3 8E 7r (k 1 )E 7r (k 2 )E a (k) 

6 (hi — k 2 — k) 



[2a t (fci)a(A; 2 )6(A;) 



-E v {k x ) + E v (k 2 ) + E a (k) 



+ <^<^ eJ^eX) + e,(^ + h - c '» ■ (2 - 15) 

By using Eq. (2.11), Eq. (2.10) can be written as 

H' = H + H' I7 (2.16) 

with 

H'j = Hf + [Hf, iS] + - [Hf, iS ] + higher order terms . (2.17) 

Since Hf, Hf, and S are all of the first order in the coupling constant g awn , all processes 
included in the second and third terms of the H'j are of the order of g^ n7r - But some of 
them are 'real processes', such as the a-exchange nil interaction and 7r-exchange ixix — > aa 
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transition , as illustrated in the upper half of Fig. 2. The other processes are 'virtual 
processes'. An example is the emission of two a mesons by a pion illustrated in the lower 
half of Fig. 2. We therefore rewrite Eq. (2.17) as 

H\ = Hf + HT + H? + Y, 0(g n an7T ) , (2.18) 

n>3 

where 

H' I p = {[Hf y iS\ + \[H?,iS\) p y (2.19) 
Hp = ([Hf,iS\ + ±[Hf,iS\f. (2.20) 

In the above two definitions, H' I P (H' I Q ) is obtained by evaluating the commutators using 
Eqs. (2.7), (2.8), and (2.15) and keep only the real (virtual) processes in the results. 

The next step is to perform a second unitary tranformation to eliminate H ®. In this 
paper we only consider the effective Hamiltonian up to second-order in the coupling constant, 
so we do not need to consider the second unitary transformation. The effective Hamiltonian 
is then obtained by dropping H'® and higher-order terms in Eq. (2.18) 

H eff = H + H P + H'f , (2.21) 

where Hf is defined by Eq. (2.7), and H' p can be calculated from Eq. (2.19) by using 
the solution Eq. (2.15) for S. Note that H'f contains a a — 7r-loop correction to the pion 
mass term. If we choose the pion mass in Hq as the physical mass, this loop correction 
should be dropped from H'f in order to avoid double counting. A similar situation will 
be encountered in the derivation of itN effective Hamiltonian. In our model we choose the 
physical masses for N and 7r in the free Hamiltonian and hence the loop corrections in the 
effective Hamiltonian are also dropped. This phenomenological procedure saves us from 
facing the complicated mass renormalization problem in solving scattering problems. 
Because of Eq. (2.14), Eqs. (2.12) and (2.13) lead to the following simple relation 

<f\iS\i>= — ^-=- </|fl?|i> • (2.22) 
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where \f > and \i > are two eigenstates of H . With the above relation, the calculation of 
H'j p Eq. (2.19) becomes 

1 



<f\H' I p \i> = E I < / \ H i \n><n \Hf \i > 

n { 



Ei — E n 



- <f\H?\n><n | iff \i > 



E„. — E 



+ <f\H?\n><n \H? \i > 1 



1 



^2.23) 



2 \_Ei — E n E n — Ef \ J 
The calculation of the matrix element of H I P therefore has a very simple rule. For a given 
choice of basis states \i > and \f >, the allowed intermediate state n is determined by the 
operator structure of O n in Eqs. (2.12) and (2.13). The denominator in Eq. (2.23) can easily 
be written down by using eigenvalues of the free Hamiltonian H . 

Evaluating Hf and H'f explicitly within the coupled nn © a sapce, Eq. (2.21) can be 
cast into the following more familiar form for im scattering 

H£ff = Hq + j~(J < — >7V7T + K-7T- (2.24) 

Here H is the free Hamiltonian operator for 7r and a mesons. The second term describes 
the a <-> mi transition with the following matrix element 

< hh | I k >= . (2.25) 

^(27rf8E 7T (k 1 )E n (k 2 )E (T (k) 

The 7T7T potential V nn is obtained by using Eq. (2.23) to calculate H'f between two tttt 
states. For \i >= \kak i2 >= \f^[ a j-. 1 a %. 2 ] |0 > an d < f\ =< kf±kf 2 \ =< ®\[ a k fl a k f2 ]\f^ the 
possible intermediate states in Eq. (2.23) are \ixita > and \ixixixixg > states. Inserting these 
intermediate states into Eq. (2.23) and earring out straightforward operator algebra, we 
obtain 

V m = V: n + Vl (2.26) 
with the following the matrix elements between two nn states 

ffLn 1111 
j aim 



< kfikfz | | k a k i2 >= 



(2^r) 3 ^2E n (k fl ) pE v {k f2 ) pE v {k tl ) j2E w (k i2 ) 

x [Di'\k a + k l2 ) + Dl-\k fl + k f2 )} , (2.27) 
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<k fl k f2 \ V^\k iX k i2 >- ^ 



where 



(2vr) 3 ^j2E n (k fl ) ^2E v (k f2 ) pE n {k iX ) pE n {k i2 ) 
x [D a (k a - k f2 ) + D a {k i2 - k f2 ) + D a (k a - k fl ) + D a (k i2 - k fl )} , (2.28) 



k z — mi 



with 



This completes the illustration of the SKO method in deriving an effective Hamiltonian 
from a model Lagrangian of relativistic quantum field theory. The extension of the method 
to consider more realistic Lagrangians is straightforward and will not be further detailed. In 
the following sections, we will simply write down the starting Lagrangians and the resulting 
effective Hamiltonians up to second order in the coupling constants for nN scattering and 
the jN — > ttN reaction. 

III. 7T-N SCATTERING 



We start with the following commonly assumed |T7[ Lagrangian for N, A, it and p fields 



L(x) = L (x) + L/(x), (3.1) 
where Lq(x) is the usual noninteracting Lagrangian, and the interaction is taken to be 

L I {x) = L 1tNN (x) + L wNA (x) + L pNN (x) + L p7r7r (x), (3.2) 



with( in the convention of Bjorkin and Drell |52[ ) 

LirNNix) = -^^^N(x)l5lvT4> N (x)d^ ■ f v (x), (3.3) 

LirNA = ^^ip A (x)fip N (x) ■ d^ix) + [h.c], (3.4) 
L pNN (x) = g P NN^N{.x)- ■ hv^pO^) - ir J! —VLLvd u (j) p \x)}il) N (x), (3.5) 

A LVTl TV 

L P 7tA x ) = 9p™(<i>* x ■ <t>p- (3-6) 
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Here T is a N — > A isospin transition operator defined by the reduced matrix element 
54fl < |||T||| >= — < |||T"I"||| >= 2. By using the standard canonical quantization, a 



Hamiltonian can be derived from the above Lagrangian except the term involving the A 
field. The difficulty of quantizing the A field is well known, as discussed, for example, in 



textbook [|55| and Ref . |56| . As part of our phenomenology, we take the simplest prescription 
by imposing the following anti-commutation relation 

{A # ,A+} = 5(p-pO, (3.7) 
where A^(A^) is the annihilation(creation) operator for a A state. This choice then leads 



[ 55] to the A propagator given later in Eq. (3.18). The alternative approaches proposed in 
Ref. [|56] will not be considered. 

Following the procedure described in section II, the next step is to decompose the result- 
ing Hamiltonian into a Hf for 'physical processes' and a Hf for 'virtual processes'. From 
the expressions Eqs. (3.3)-(3.6), it is clear that the real processes in this case are A «-> irN 
and p «-> 7T7T transitions which can take place in free space(because > + m n and 
m p > 2m n ). The virtual processes are N <-> ttN, N <-> pN ,N <-> 7rA, and tt <-> Tip tran- 
sitions. These virtual processes can be eliminated by introducing a unitary transformation 
operator S which can be determined by using the similar method in obtaining the solution 
Eqs. (2.13)-(2.14). Here, we of course encounter a much more involved task to account for 
the Dirac spin structure, isospin, and also the anti-particle components of N and A. To see 
the main steps, we present in Appendix A an explicit derivation of the potential due to the 
LnNN term Eq. (3.3). 

For practical applications, it is sufficient to present our results in the coupled ttN © 
A subspace in which the ttN scattering problem will be solved. The resulting effective 
Hamiltonian then takes the following form 

^eff — Hq + Ta^ttAT + v wN , (3.8) 

where H is the free Hamiltonian for tt, N and A. Note that T^^n (Figs. 3a and 3b) is 
the only vertex interaction in the constructed effective Hamiltonian. Our model is therefore 
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distinctively different from the previous meson-exchange irN models |T^-[T4|,[T9[ which all 
involve a bare nucleon state N and aiVo<-> ttN vertex. 
The ttN potential v n N in Eq. (3.8) is found to be 

VnN = v Nd +v Ne +v p + va d +v Ae , (3.9) 

where v^ D is the direct nucleon pole term (Fig. 3c), v^ E the nucleon-exchange term (Fig. 
3d), v p the p— exchange term (Fig. 3e), v& D the interaction due to the anti-A component 
of the A propagation (Fig. 3f), and va e the A-exchange term (Fig. 3g). To simplify the 
presentation, we will only give the matrix element of v^n in the ttN center of mass frame. 
The initial and final four-momenta k^, k'^ for pions and p^,p'^ for the nucleons in Fig. 3 are 
therefore defined as 

k» = (E„(k),k), 
^ = {E N {k),-k), 
k' 11 = {E n {k'), k), 

p'» = (E N (k'),-k). (3.10) 

In terms of these variables, the matrix element of each term of Eq. (3.9) between two ttN 
states can be written as 



< k % , m s m T \v a \ki, m s , m T > 



( 27r ) 3 JlEjtt)\ E N(k') pE„(k)\ E N(k) ' 
x u_k, tm , tm , T I a (k'i'i ki)u_ %ms>mT , (3.11) 

where u^ m3i m T is the Dirac spinior, m s and m T are the nucleon spin and isospin quantum 
numbers, i and i' are the pion isospin components. The interaction mechanisms are con- 
tained in the functions I a (k'i', ki). After performing lengthy derivations, we find that these 
functions can be written in the following concise forms 

I ND (k'i',ki) = (^)W ^[S N (p + k) + S N (p' + k')]ra 5 Jf, (3.12) 

I NE (k'i', ki) = (^) 2 r i7 5 #hs N (p - k') + S N (p> - k)]W Jt>, , (3.13) 
m-x 2 
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I p {k'i', ki) 



Ja d (^V, ki) 



lAE(k'i', ki) 



igpNNg P TTTT f r k„ 

: ^u'kTki [7m - 



i^(p-pT}Df(p- P ')(k + k') x ) 

Z TYl f 



+[(p-p')~(k'-k)]}, 



1 



TK-[S^ip + k)+S^{p' + k') 



S^ip + k) 



S^ u (p' + k')}T t k u 



■fir 



m 



—" 2 T}k^[S%{p - k') + - k)\T v k' v 



(3.14) 

(3.15) 
(3.16) 



The propagators in the above equations are defined as 
1 



S N (p) 
5T(p) 

£>r(p) 



1 



3(^- m A ) 



2(-<T + ^f) + 



P I 



m. 



g^ u — p^p u /m 2 
p2 _ m 2 



(3.17) 
(3.18) 
(3.19) 



In Eq. (3.15), we also have introduced a propagator 



S 



(P) 



v v 



Ea(p) Po ~ E A (p) 
where uj£ is the Rarita-Schwinger spinor(as explicitly defined in Ref. Jl 
frame, this propagator reduces to the following simple form 



S( +)il {p) 



1 + 7 C 



(3<y*- 



1 



(3.20) 
In the A rest 

(3.21) 



p —y 6 po — m& 

for i, j = 1,2, 3. The other elements involving time components vanish in this special frame; 
g(+)/i0 _ g(+)ov _ q ^Yie appearance of this propagator in Eq. (3.15) is to remove the 
ttN — > A — > ttN mechanism which can be generated by the vertex interaction Ta^ttN of the 
effective Hamiltonian Eq. (3.8). This comes about naturally in our derivations. 

We note that the above expressions are remarkbly similar to those derived from us- 
ing Feynman rules. The only differences are in the propagators of the intermediate par- 
ticles. These propagators are evaluated by using the momenta of the external particles 
which are restricted on their mass shell, as defined in Eq. (3.10). For the off-energy-shell 
dynamics(i?Ar(fc) + E n ( k) ^ E N {k') + E^ik')) , these propagators can have two possible forms, 
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depending on which set of external momenta is used. The propagators in Eqs. (3.12)-(3.16) 
are the average of these two possible forms of propagators. More details can be seen in 
Appendix A where the derivation of vn = vn d + vn e is given explicitly. 

In the irN center of mass frame, the A in the vertex interaction Ta^ttAt is at rest. In 
this particular frame, the Rarita-Schwinger spinors reduces to a simple form such that the 
matrix element of the vertex interaction Ta^ttN takes the following familiar form 

f-rrNA i 1 



< A I^Wn-AT \ki > 



EN{k) + mN S.k Tl . (3.22) 



'{2-Kf y/2Ejjz) \ 2E N {k) 
Here S is a. N —>■ A transition spin operator. It is defined by the same reduced matrix 
element as the transition isospin operator T. 

Because of the absence of a ttN «-> N vertex in the effective Hamiltonian Eq. (3.8), 
it is straightforward to derive the ttN scattering equations in the coupled ttN © A space. 
The derivation procedure is similar to that given in Ref. |23] for the more complicated 7rNN 



problem. The essential idea is to apply the standard projection operator technique of nuclear 
reaction theory ||. The resulting scattering amplitude can be cast into the following form 

T^n^E) = t wN (E) + T a ^ 7tN (E)Ga(E)T 7tN _^ a (E). (3.23) 

The first term is the nonresonant amplitude determined only by the irN potential 

t nN (E) = v nN + v 7rN G 7rN (E)t nN (E), (3.24) 

with 

E — E N (k) — E n (k) + xe 
where P n ^ is the projection operator for the nN subspace. The second term of Eq. (3.23) is 
the resonant term determined by the dressed A propagator and the dressed vertex functions. 
They are defined by 



G A (E) = G° A (E) + G° A (E) E a (E) G A (E) 

- Pa 
E — m A — T, A (E) ' 
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(3.26) 



with 
and 

TwN-^a(E) = r 7rA r__ } A(l + C 7r Ar(-E)t 7r Ar(-E)), (3.28) 

(E) = (1 + GW(£)^v(£))r A ^v, (3.29) 

where Pa is the projection operator for the A state, and the A self-energy is defined by 

(25) . (3.30) 

Equations (3.23)-(3.30) are illustrated in Fig. 4. These equations are solved in partial- 
wave representation. To find the solution for the integral equation (3.24), it is necessary to 
regularize the 7rN potential by introducing a form factor for each vertex in Eqs. (3.12)-(3.16). 
In this work, we choose 

[Ind + lN E }(k'i', ki) -> [I ND + I NE \{k'i'i ki)F^ NN {k')F^ NN {k) , (3.31) 
I p (k>i\ki)^I p (k7,^F pNN (k-k')F pnn (k-k'), (3.32) 

(3.33) 



with 



F, NN (k) = ( ^ 2 AlNN ' 2 ) 2 , (3.34) 

F pNN (k - k') = F p7Tn (k - %') , (3.35) 

= ( \ p - f. (3.36) 

A2 + (A;-A;') 2 

For irNA vertex with an external pion momentum k, we choose 

F,NA(k) = ( .f^r/ - (3-37) 

We have also tried other parameterizations of form factors, but they do not give better fits 
to the tcN scattering phase shifts. 
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IV. PION PHOTOPRODUCTION 



To proceed, we need to first extend the Lagrangian Eq. (3.1) to include uj meson coupling 
which is known PJT7f to play an important role in pion photoproduction. We choose the 
following rather conventional form (with ~ 0) 

L^nn = 9wNNi/>N(x)bM(x) - —^—a liV d v (j)^{x)\^ N {x). (4.1) 



Following the approach of Ref. [TjJ, the pion photoproduction mechanisms are defined by 



the hadronic Lagrangians defined by Eq. (3.1) and (4.1) and the following electromagnetic 
interaction Lagrangians 

with 

L-yNN = -e$ N {x)[eA(x) - -^—a^A^x))]^^), (4.3) 

L-ynNN = - C ^ NN ip N (x)j 5 A(x) (f x $ 7r (x)) 3 ip N (x), (4.4) 

L 77r7r = e(<9 M 7r (x) x ^(x^^A^x), (4.5) 

L p7r7 = i^Le a p y s(d a Af(x))Mx) ■ (d^%x)), (4.6) 

L_ 7 = 9 ^e aMS {d^A\x))4>l{x){d^i{^ (4.7) 
m n 

L-yNA = iej%T 3 T lu/ A u (x)il> N (x) + [h.c.]. (4.8) 
Here e = (1 + r 3 )/2, £ = (k p + « n )/2 + (/%, — n n )r 3 /2. The 7iVA coupling in Eq. (4.8) is 

T» v = {G M K™ + G E Kl), (4.9) 
as defined in Eqs. (2.10b,c) of Ref. [I7|. Its matrix element between an iV with momentum 



p and a A with momentum p& can be written explicitly as 

< A( PA )\T, u \N(p) >= (G M - G E )[- 3 — - mA +mN 6 M ^PV] 

(wa — m N y — q 2 2m N 

+<W[t ? V2 ^7 V2 m * + mN (4.10) 

((m A + mjv) 2 - q 2 ){(m A - m N ) z - <r) m-Ar 1 
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with P = (p + Pa)/2 and p& = P + <?• 

By appling the usual canonical quantization procedure, we can obtain from the above 
Lagrangians an electromagnetic interaction Hamiltonian H em . In this work, we will treat 
the electromagnetic field as an external classical field, and hence the electromagnetic inter- 
action H em can be neglected in constructing the unitary transformation operator S. The 
effective Hamiltonian for describing pion photoproduction is therefore a simple extension of 
the effective Hamiltonian of the form of Eq. (2.21) 

U v U i Tjem 

H eff -> ti e ff + H eff 

= H + Hf + H'f + Hlf f (4. 11) 

with 

H^ f = H em + [H em ,iS] (4.12) 

By evaluating Hf^f in the coupled A © nN © 7iV subspace, we obtain an extension of Eq. 
(3.8) 

Heff H]ff — H + Ta^ttN + VwN + Ta^ 7 at + V ny , (4-13) 

where T^ nN and v^n have already been given in Eqs. (3.1 1)-(3.22). The resonant and 
nonresonant electromagnetic interactions are, respectively, described by Ta^tv and v ni , 
and are illustrated in Fig. 5. We again omit the details of the derivation of these two terms, 
and simply present our results in the center of mass frame. The momenta variables for 
the photon, p^ for the initial nucleon, k^ for the pion, and p' M for the final nucleon in Fig. 5 
are therefore 

<f = (q,q), 

p» = (E N (q),-q), 

p> i = (E N (k),-k). (4.14) 
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In terms of these variables, the expression Eq. (4.10) becomes very simple. The matrix 
element of the T^ lN vertex can then be expressed in terms of the spin and isospin iV — > A 
transition operators S and T introduced in section III. At the resonance energy, m A = 
E N (q) + q, Eq. (4.10) leads to 



< A|r 7A ^ A |gA > 



E N (q) + m N 1 3m A 



(2tt) 3 / 2 \| 2E N (q) ^/2q~ 2m N (m A + m N ) 
T 3 [iG M S x q. e x + G E (S ■ e x a ■ q + S ■ qa ■ e A )], (4.15) 



where e\ is the photon polarization vector. The matrix element of the nonresonant interac- 



tion Vn-y can be written as 



r- i /. i i 1 / rn N I m N 

< ki,m s m T \v ni \qX,m s m T > = 



x S-^EW&^hM^- (4.16) 

a 

The nonresonant pion photoproduction mechanisms are contained in I 1 ^ for the direct nu- 
cleon terms (Figs. 5c,5d,5e) , for the pion pole term (Fig. 5f), ITJ^ for the vector meson 
exchange (Fig. 5g), and /ad,A£ for the direct and exchange A terms (Figs. 5h,5i). Explicitly, 
we have 

qX) = e -^[irrf $S N {p' + k)(e /(A) + i— ^— <7 /1 „e(A) / V) 
777^ 2m w 

+,(e /(A) + i_^_< V e(A)'V)Siv(p - A;)r l7 5 ft - e ij3 r n b /(A)] (4.17) 

ir(U,q\) = -^Le ij3 rrf(p"- Mk + p-p')-e(X)D n (p-p') (4.18) 
I?(lk, q\) = ^^[7, - ^-a„(p - p'T] x D?{p - p')e a ^(X)q^p - p'V 

I?(ki, qX) = 9uNn9u ^ t 3 [^ - ^-a m (p - p'Y>] x D^(p - p')e a ^e a (X)q^p - p'V (4.19) 

Ill(ki, qX) = e -^TlT\ s e\X)S^{p - k)T,k v (4.20) 

lH(ki, qX) = - e J^T}k,[S%{p' + k) - S^(p' + k)]T 3 T uS e 5 (X) (4.21) 

HI:. 

Here we observe again that the above expressions are very similar to the results derived by 
using Feynman rules. However, they have an important feature that the time components 
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of the momenta in the propagators and strong interaction vetices are evaluated by using 
the external momenta of the final nN state. This is the consequency of appling the unitary 
transformation method defined in Eq. (4.12). In addition to including nonresonant A terms 
(Figs. 5h and 5i), this is an another feature which makes our model different from the model 
developed in Ref. [ITT . 



It is straightforward to derive from the effective Hamiltonian Eq. (4.13) the t-matrix of 
pion photoproduction 

T yn (E) = t J7r (E) + r A ^ N (E)G A (E)T lN -+A(E) , (4.22) 

where the nonresonant amplitude is defined by 

t 7n (E) = Vrps + t nN {E)G^ N {E)v y7r . (4.23) 

The dressed <-> A is defined by 

^^n^a(E) = r 7 7v^A + F 7tN ^a(E)G 7tN (E)v 17T . (4.24) 

In the above equations, G A , T A ^ nN , G^n an d t nN have been defined in section III. The 
standard partial-wave decomposition is used to obtain the multipole amplitudes from T 77r 
for the 7iV — > ttN reaction, and from T A ^ n N for the dressed A «-> ttN vertex. Eqs.(4.22)- 
(4.24) are illustrated in Fig.6. 

The K matrix formulation of the 7iV — > nN reaction is often used in the analysis of 
data. Within our formulation, this can be obtained by replacing the irN free Green function 
G, N Eq. (3.25) by 

GME)^GUE)=P E _ EN ^_ K(ky (4.25) 

where P means taking the principal-value part of the propagator. If this replacement is used 
in the calculations of Eqs. (3.23)-(3.30), all scattering quantities will be real numbers. These 
i^-matrix quantities are defined by exactly the same Eqs. (3.23)-(3.30) with the following 
changes 
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G-nN ~~ > G^ N , 
T n N — > KnN , 
tnN * k n N , 

with 

£ A (£) = r wJV _ A G^(S)rl^(S) . (4.27) 

Note that the A self-energy £ A is now a real number, and the propagator G A has a pole at 
£ = M* = m A + E A (M R ). 

The corresponding K-marix for pion photoproduction can be obtained from Eqs. (4.22)- 
(4.24) by the same replacement Eq. (4.25) 

K^{E) = k Jn (E) + T k A ^ N (E)G p A (E)T k N _ A (E) , (4.28) 

with 

k 17 r = *V + Kn(E)G p n (E)v^ , (4.29) 

rV» A (£) = r 7iV ^ A + r k N _^ A (E)G P N (E)v, n . (4.30) 

For the on-shell matrix elements(E = E N (k ) + E w (k ) = q + E N (q)), it is straigtforward to 
find the following relation in each partial wave 

T-y W (k , q) = (1 - mpT„ N (ko, k ))K^(k , q) . (4.31) 

For investigating the hadron structure, we are interested in the 7iV <-> A vertex. As 
seen in Eqs. (4.24) and (4.30), the dressed vertices in t-matrix and in K-matrix are different. 
In the t-matrix formulation r 7 jv )A Eq. (4.24) is a complex quantity, while in the K-matrix 
formulation T k N A , Eq. (4.30), is a real function. Consequently, we need to be careful about 
the meaning of the E2/M1 ratio of the dressed 'jN <-> A vertex. The clearest definition 
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seems to be in the K-matrix formulation because as the energy approaches the resonance 
position M R = m A + E A (M R ), Eq. (4.28) is reduced to 

with 

A = T A ^ N (E)T k N ^ A (E) (4.33) 
B = k^{E) (4.34) 

The separable form of the residue A of the K-matrix leads to an interesting result that 
the ratio between the El and Ml multipole amplitudes of the dressed 7AA vertex can 
be directly calculated from the residues of the corresponding multipole amplitudes of the 
7iV — > nN reaction. The reason is that both amplitudes have the same strong interaction 
dressed vertex in the P 33 channel, and hence the ratio between the residues does not depend 
on it. Explecitly, we have 



Rem ~ l A(M^) hN ^ N 



i(M a 

ri^v(P 3 3)r^ A (£ 1+ ) 

fi_^(P33)f^ A (M 1+ ) 

T 7N ^ A (E 1+ ) 



(4.35) 



T lN ^ A (M 1+ ) 

The above relation is the basis of the model-independent analysis of Ref. |J. We will discuss 
this issue in the next section. 



V. RESULTS AND DISCUSSIONS 

Our first task is to determine the parameters of the effective ttN Hamiltonian derived in 

f' 2 

section II. Apart from the known irNN coupling constant, - n ff L = 0.08, the model has seven 
parameters: the couplig constants {gl = Q pN n Q pnn ? k p , fnNA) of Eqs. (3.12)-(3.16), the cutoff 
parameters (A^nn, ^na, A p ) of the form factors Eqs. (3.34)-(3.37), and m A of the A bare 
mass. These parameters are determined by fitting the irN phase shifts. Without including 

24 



inelastic channels, our scattering equations Eqs. (3.23)-(3.30) are valid rigorously only in the 
energy region where the nN scattering is purely elastic. We therefore first take a conservative 
approach to only fit the data in the energy region below Tl = 250 MeV pion laboratory 
energy This model, called Model-L, is sufficient for investigating pion photoproduction up 
to 400 MeV photon laboratory energy. Our results are displayed in Fig. 7. We see that 
within the uncertainties of the phase shift data f38H4"0| the model can give a good account of 



all s- and p- partial waves except the P13 channel at Tl > 120 MeV. We have found that this 
difficulty can not be removed by tring various form factors other than those given in Eqs. 
(3.34)-(3.37), and following the previous works fl2j|l3l to include the exchange of a fictituous 



scalar a meson. To see the origin of this problem, we show in Fig. 8 the contributions from 
each mechanism of Fig. 3 to the on-shell matrix elements of the ttN potential. Clearly, 
the fit to the phase shift data involves delicate cancellations between different mechanisms. 
It is possible to imporve the fit to P13 by weakening the p-exchange or the A-exchange. 
But this change will destroy the good fits to all other partial waves. Fortunately, the ttN 
scattering effect due to the P13 channel is weak in determining the pion photoproduction 
cross sections. We therefore will not pursue the solution of this problem here. Perhaps this 
can be solved only when the p-exchange is replaced by the two-pion-exchange considered 



in Ref. ||14j| . To be consistent, the coupling with two-pion channels, such as irA and pN, 
must also be included. These two possible improvements can be achieved by extending the 
unitary transformation method introduced in Sections II- IV to second-order in the coupling 
constants. 

Let us now examine in more detail the P33 channel which is most relevant to our later 
investigation of the A excitation in pion photoproduction. As seen in Eq. (3.23), the resonant 
part of the P 33 amplitude is determined by the dressed propagator of Eq. (3.26) and 
the dressed vertex f a<->kN of Eqs. (3.28) and (3.29). Clearly, the A resonance peak of ttN 
scattering can be obtained only when the model can generate a A self-energy such that 
the real part of (E — — Sa(P)) — > as the ttN invariant mass W approaches the 
resonance enery W = Mr = 1236 MeV. Our model has this desired property, as illustrated 
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in the upper half of Fig. 9. Another important feature in the P33 channel is that the itN 
potential generates the dressed A <-> ttN vertex, as defined in Eqs. (3.28) and (3.29). We 
have found that this renormalization effect modifies greatly the A <-» nN form factor in the 
low momentum region. To see this, we cast the bare vertex (Eq. (3.22) including the form 
factor F nN &(k) defined by Eq. (3.37)) and the dressed vertex (Eq. (3.28)) into the following 
forms 

(A\T A ^ N \H) = - fwNA % —^=F bare (k)S - kTi. (5.1) 

^(2^)3 v / 2^ 

(A|f A^jvl^*) = -hlE^L_^ = J; — F dressed (k)S ■ kT^ (5.2) 



with the normalization \Fd reS sed(Q)\ = Fb are (0) = 1. We find that the dressed coupling 
constant, f n NA, is 1.3 of the bare coupling constant f n NA- The dressed form factor F n NA(k) 
falls off faster than the bare form factor F n ^^(k) in momentum space, as seen in the lower 
half of Fig. 9. This means that the nonresonant tcN interaction has extended the A excitation 
region to a larger distance in coordinate space. 

A significant difference between our approach and the previous ttN models |I12|-|Il|,|r9|] is 



in the treatment of Pn channel. By employing the unitary transformation, the nN <-> iV 
vertex does not appear in the effective Hamiltonian Eq. (3.8) and hence our formulation of 
ttN scattering is straightforward. It does not require the nucleon mass renormalization. It 
is natual to ask whether our approach possesses the well-established nucleon-pole dynamics. 
This question can be answered by examining Fig. 10 in which the irN phase shifts and the 
scattering t-matrix elements calculated from the nucleon pole term vn d only (Fig. 3c) and 
the full potential are compared. We see in the the upper half of Fig. 10 that ttN phase shifts 
due to the nucleon-pole term (dotted curve) are repulsive as expected. The fit to the Pn data 
is due to a delicate cancellation between the repulsive nucleon-pole term and the attraction 
coming mainly from p— and A— exchange terms(see the Pn case in Fig. 8). In the lower half 
of Fig. 10, we see that as the energy approaches the threshold, W = m n + mjy, the nucleon 
pole term(dotted curve) apparently dominates the interaction. If we analytically continue 
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to the nucleon pole position, k = ik^im^ = E^iikN) + E n (ik]y)), the scattering amplitude 
will be determined by the nucleon pole term, i.e. t(k , k ,W — > m^) ~ f jv d ~ 1/(W — m iv)- 
The parameters of the constructed model are listed in the first row (Model-L) of Table 
I. The calculated scattering lengths are presented in the first column of Table II. They are 
all in good agreement with the data. If we assume the universality of p coupling, we then 
have g P NN = g P mr = 6.2 which is close to that determined in Refs. |12|,|13| . The fit is also 
sensitive to the p tensor coupling constant n p . Our value is close to that of Ref. jl2|, but is 
much smaller than 6.6 used in Ref. |13 . 



We now turn to presenting our results of pion photoproduction. With the ttNN, it AN 
and pNN vertices defined by the parameters given in Table 1, the considered pion photo- 
production mechanisms (Fig. 5) still have unknown parameters associated with the vector 
meson-exchange and the 7iV <-> A vertex. Following the previous approach |17|| , we assume 
that the photon-meson coupling constants g pnl and <7 W7r7 can be determined from the partial 
decay widths listed by the Particle Data Group For the uj meson, we further assume 
that the tensor coupling k^nn — and the ujNN form factor is identical to the pNN form 
factor given in Table 1. The coupling constant g^NN is not well determined in the literature. 
We will treat it as a free parameter, although the quark model value g W NN — (^9pNn)/^ 
seems to be a reasonalbe guess. Thus, our investigation of pion photoproduction has only 
three adjustable parameters: Gm and Ge of the bare A <-> vertex, and the coupling 
constant g U NN of the uj exchange. We have, however, some ideas about the ranges of these 
parameters. If we assume that bare vertex interaction Ta^ 7 v can be identified with the 



constituent quark model |^6| , ^8| , ^7| , ^9| , |50| , then \Ge/Gm\ ~ since the one-gluon-exchange 
interaction gives negligible -D-state components in iV and A. We also expect that the uj 
coupling should be close to the quark model prediction, g u NN = ^g P NN/^ ~ 9, if the p 
coupling from our irN model (Table I) is used. It is therefore reasonable to only consider 
the region g^NN < 15 and \Ge/Gm\ < 0.1. 

Since the tu-exchange mechanism (Fig. 5g) does not produce charged pions directly (only 
through Ti N charge exchange), the ranges of g u NN, Gm and Ge can be most sensitively 
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determined by considering the data of ir° photoproduction. In the considered region that 
\Ge/Gm\ < 0.1 and g U NN — 15, we have found that the magnitudes of the 7r° differential 
cross sections depend mainly on Gm- The values of g W NN and Ge can be narrowed down 



by considering spin observables. In this work, we make use of the recent LEGS [[ITJ data 
of photon-asymmetry ratios R 1 = dan/ da ± of 7^ — > 7i°p reaction. We have found that the 
slope of Ry(E) at a fixed pion angle is sensitive to the value of g^NN ■ This is illustrated in 
the upper half of Fig. 11 for the case of Gm = 1-85 and Ge = 0.025. A smaller g^NN yields 
a steeper slope. The data clearly favor g^NN ~ 10.5. In the lower half of Fig. 11, we see 
that the magnitude, not the slope, of _R 7 is significantly changed by varing the value of Ge 
from -0.1 to +0.1. The data are consistent with -0.025 <G E < 0.025, while G E = +0.025 
seems to give a better fit. Results similar to Fig. 11 can be obtained by using a higher value 
Gm = 1-95. In this case, a smaller value of g^NN = 7 is needed to maintain the same fit to the 
magnitude of the differential cross section as well as the slope of the R 1 . But the best value 
of Ge to reproduce the magnitude of R 1 is —0.025 instead of +0.025 for the g^NN — 10-5 
case. In fact, we have observed a strong correlation between the allowed values of and 
g^NN- In all cases, the allowed value of Ge is consistent with —0.025 < Ge < +0.025. 
Therefore, the acceptable values of (Gm, <?u>aw) are on the curve between the Ge = —0.025 
and Ge = +0.025 lines in Fig. 12. To determine the precise value of Rem = q^, which 
measures the deformation of the A, we clearly need to pin down the u meson coupling 
constant g^NN- 

The predictions from using the parameters lying on the curve between Ge = +0.025 and 
Ge = —0.025 lines of Fig. 12 are also in good agreement with the 7r° data at other angles and 
the 7r + data. These are shown in Fig. 13 for tc° production and Fig. 14 for tt + production. 
These results are obtained from using the parameters defined by the interaction points of the 
G E = ±0.025 lines in Fig. 12: (g uNN , G M , G E ) = (10.5, 1.85, +0.025), and (7., 1.95, -0.025). 
Both sets of parameters yield equally good agreements with the ir° data (Fig. 13). For 7r + 
production, the predictions are in good agreements with the data of the photon-asymmetry 
ratios _R 7 , but underestimate the differential cross sections by about 10 percent at most 
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energies. Since the u exchange has a small contribution to the 7r + production (only through 
charge-exchange nN final state interaction), the only way to resolve this difficulty within 
our model is to increase the value of Gm- But this will lead to an overestimate of the 7r° 
cross section from Bonn. 

Although the difficulty in reproducing the 7r + data in Fig. 14b could be an indication of 
the deficiency of our model, the possibility of a larger 7r° cross section has been suggested by 
three ir° data at 9 = 120 from Ref. |44| (Fig. 13b). To fit these three data points, we need 
to increase Gm from 1.85 to 2.0 for the case of g U NN — 10.5 and from 1.95 to 2.1 for the 
case of g W NN = 7. The results from these two changes in Gm are, respectively, the solid and 
dotted curves in Figs. 15 and 16. The predicted photon-asymmetry ratios (Figs. 15a and 
16a) are still in good agreement with the data. The agreement with the tt + data (Fig. 16b) 
is clearly improved. But the calculated ir° differential cross sections (Fig. 15b) overestimate 
the Bonn data |42]j43[1 by about 15 percent. Clearly, the disagreement between the n° data 
at 9 — 120° (Fig. 15b) from Refs. f42] , |43"| and |0| must be resolved by new measurements. 



To further reveal the dynamical content of our model, we compare in Fig. 17 our 
predictions of angular distributions with the data compiled in Ref. |4l| for jp — > ir°p 
(Fig. 17a), 7p — > 7r + n (Fig. 17b), and — > 7t~p (Fig. 17c) reactions. All results cal- 
culated from using the parameters lying on the curve between the Ge = +0.025 and 
Ge = —0.025 lines in Fig. 12 are very close to the solid curves which are from using 
(9wNN,Gm,Ge) — (10.5, 1.85,0.025). Again, we see that the charged pion production cross 
sections are underestimated. If a larger Gm = 2.0 is used in this calculation, we obtain the 
dotted curves which are in a better agreement with the charged pion data, but overestimate 
the 7T° data by about 15 percent at resonance peaks. In all cases, the theoretical predictions 
underestimate the data at 380 MeV and higher energies. This is expected, since the con- 
structed model does not include inelastic channels which should start to play a significant 
role at energies above about 350 MeV. For example, the inelastic production mechanism 
7iV — > 7i A — > ttN should exist since it is known that the ttN scattering at this higher energy 
can be described only when the coupling with the ttA channel is included. To investigate 
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this effect, it is necessary to extend the derivation of effective Hamiltonians presented in 
sections III and IV to include the 7rA as well as other two-pion states. 

We now focus on the theoretical interpretations of the A <-> 7iV vertex. The values of Gm 
and Ge determined above characterize the bare A <-> 'jN vertex which can only be identified 
with hadron models with no coupling with the ttN or other hadronic reaction channels. One 
possible interpretation is to compare the determined Gm and Ge with the predictions of the 
most well-developed constituent quark model |^6j^8| , ^7| , ^ , p0i . To explore this possibility, 



it is necessary to first discuss the quantities in our model which can be compared with the 
results from empirical amplitude analyses ||[7|,|57| ■ For investigating the A mechanism, we 
need to only consider the 7^ — > ttN multipole amplitudes M]+ and E\+ with a P33 final 
ttN state, and the dressed vertex function f 7 jv,a- These can be computed from Eqs. (4.22)- 
(4.24) or Eqs. (4.28)-(4.30) by performing the standard partial-wave decomposition (see, 
for example, the appendix of Ref. [fL7|). We will discuss these quantities using the results 
calculated from setting (g^NN, Gm, Ge) = (10.5, 1.85, 0.025) (solid curves in Figs. 13, 14 and 
17). 

The predicted amplitudes M\+ and are compared in Fig. 18 with the results from 
the empirical amplitude analyses Ref. |7|,|57|. We see in the upper part of Fig. 18 that the 



predicted M\+ amplitudes are in good agreement with empirical values. In the lower half, we 
show that both the E x + amplitudes calculated from using Ge = +0.025 (solid curves) and 
Ge = —0.025 (dotted curves) are within the uncertainties of the amplitude analyses. This 
is consistent with our analysis using LEGS data, as seen in the lower part of Fig. 11. The 
uncertainties of the empirical values of the E\+ amplitude are due to the lack of complete 
data of spin observables. More experimental efforts are clearly needed to pin down the value 
of Ge which is needed to test models of hadron structure. 

The dressed A <-> 'jN vertex, defined by Eq. (4.24), is a complex number. By making 
the usual partial-wave decomposition, its magnetic Ml and electric E2 components can be 
written as T(a) = \T(a)\e l ^ a > with a = M 1 +,£' 1 +. The predicted dressed vertex functions 
r(a) are the solid curves in Fig. 19. We see that their magnitudes |T(a)| are very different 
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from the corresponding values(dotted curves) of the bare A <-» ^yN vertex. The differences 
are due to the very large contribution of the nonresonant mechanism described by the second 
term of Eq. (4.24). Our results indicate that an accurate reaction theory calculation of the 
nonresonant pion photoproduction mechanisms is needed to determine the bare A <-> jN 
vertex from the pion photoproduction data. This requires a dynamical treatment of the 
nonresonant pion photoproduction mechanisms, as we have done in this work. Within the 
meson-exchange formulation presented in this work, the determined Gm arid Ge of the bare 
A <-> 7iV vertex can be compared with the predictions from a hadron model which does not 
include the coupling with the irN "reaction" channel (both pion and nucleon are on their 
mass-shell). 

We now turn to investigate the X-matrix method which has been the basis of the em- 
pirical amplitude analyses of Refs. @,[7|]. In Ref. ||, it was shown that if the background 
term is assumed to be a slowly varing function of energy, the ratio Rem between the E\+ 
and Mi+ of the 'jN — ► A transition at the resonant energy W = Mr can then be extracted 
"model-independently" from the data of the ^yN — > ttN reaction. The only limitation is the 
accuracy of the employed ttN amplitudes and the Mj+ and Ei+ multipole amplitudes of the 
7iV — > ttN reaction. By using Eqs. (4.32)-(4.35), we can examine whether the i^-matrix 
method of Ref. || is consistent with our dynamical model. In Fig. 20, we display our pre- 
dictions of the energy-dependence of the total K- matrix (solid curve), the contribution from 
the resonant term (dotted curve) which has a pole at W = 1236 MeV, and the nonresonant 
contribution B (dashed curve). Clearly, the energy- dependence of the nonresonant term B 
is rather weak. The assumption made in the empirical analysis of Ref. || is fairly consistent 
with our dynamical model. 

By using Eq. (4.33), we can calculate the residue A of the i^-matrix from the dressed 
vertex r^ iV _ >A defined by Eq. (4.30). The results (solid curve) for the M\+ transitions are 
compared with that calculated from using the bare vertex T 7 at^a in the the upper half 
of Fig. 21. Similar to the results in Fig. 19 in the t-matrix formulation, we see the large 
nonresonant mechanisms in dressing the 'yN — > A vertex. The corresponding E2/M1 ratios 

31 



Rem are compared in the lower half of Fig. 21. The nonresonant mechanisms change the 
ratio by a factor of about 2 at resonance energy W = 1236 MeV. 

In Table IV, we list the predicted E x + and M x + amplitudes of the A 'jN vertex 
evaluated at the resonance energy W = 1236 MeV. The parameters {g^NN, Gm, Ge) = 
(10.5, 1.85, 0025) and (7., 1.95, -0.025) from the fits to the data (Figs. 13 and 14) are used 
in this calculations. We see that our average value Rem = (-1-8 ± 0.9)% is not too different 
from the average value (-1.07 ± 0.37)% of the empirical analysis f|. Since the assumption 
made in Ref. || is consistent with our model as discussed above, the difference perhaps 
mainly comes from the experimental uncertainties of the multipole amplitudes employed in 
the analysis of Ref. ||. The differences between our predicted multipole amplitudes and the 
empirical values shown in Fig. 18 could also be responsible to this discrepancy. To compare 
our results with the values listed by the Particle Data Group (PDG) ||45|| , we calculate the 
helicity amplitudes by 

A 3/2 = ^-[E 1+ - M 1+ ] , 
A 1/2 = --[3£ 1+ +M 1+ ]. 

The results at the resonance energy W = 1236 MeV are listed in Table V. The predictions 
from two constituent quark models ||48"| , f47H are also listed for comparision. We notice that 



our bare values are close to the constituent quark model predictions ||7|,|48| , and the dressed 



values are close to the values of PDG [[|IJ. This suggest that our bare vertex can be identified 
with the constituent quark model. The long-standing discrepancy between the constituent 
quark model predictions and the PDG values is due to the nonresonant meson-exchange 
production mechanisms which must be calculated from a dynamical approach. Similar 
considerations must be taken in comparing the PDG values with the predictions of higher 
mass N* resonances from hadron models. 

The results we have presented so far are based on the ttN model determined in a fit 
to the ttN phase shifts only up to 250 MeV. It is interesting to see the extent to which 
this model can be extended to a higher energy region where the inelastic processes are still 
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not dominant. More importantly, we would like to examine whether the extended model 
can yield significantly different nN off-shell dynamics which perhaps can help resolve the 
difficulty in reproducing the magnitudes of ir + cross section (see Fig. 14b). To explore these 
possibilities, a ttN model is constructed by fitting the phase shifts data up to 400 MeV. This 
model is called Model-H to distingush it from Model-L from the fit up to only 250 MeV. The 
resulting parameters are also listed in Table 1. The phase shifts calculated from these two 
models are compared in Fig. 22. We see that Model-H (dotted curves) clearly gives a much 
better fit to the data in the entire considered energy region. But it is not as accurate as 
Model-L in describing the the crucial P33 channel in the low energy region. To accurately fit 
the P33 in the entire energy region and to resolve the difficulty in the P± 3 channel, additional 
mechanisms may be needed. 

The — > nN results calculated from using the Model-H and Model-L are compared 
in Fig. 23. The photon-asymmetry ratios (Figs. 23a, 23c) are equally well described by both 
models. They yield, however, significant differences in describing the differential cross sec- 
tions. In Fig. 23b, we see that Model-H gives a much better description of the n° differential 
cross sections in the high energy region. But it slightly overestimates the cross sections at 
low energies. The ir + differential cross sections are better described by Model-H, as seen in 
Fig. 23d. But the difficulty in reproducing the magnitude is not removed entirely. 

The results in Fig. 23 suggest that our predictions do depend to some extent on the 
accuracy of the constructed irN model in describing the irN phase shifts. A natural next step 
is to extend the present model to include the inelastic channels to obtain an accurate fit up 
to 400 MeV. This extension then will introduce inelastic pion photoproduction mechanisms, 
such as the 7-/V — > n A — > nN process, which may be needed to resolve the difficulty in 
getting an accurate description of both the 7r° and n + processes. Such a coupled-channel 
approach must also include the effect due to the excitations of higher mass iV* nucleon 
resonances. This must be pursued in order to make progress in using the forthcoming data 
from CEBAF to test hadron models. 
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VI. CONCLUSIONS AND FUTURE STUDIES 



We have applied the unitary transformation method first proposed in Ref. [22| to derive 
from a model Lagrangian with N, A, 7r and p, u, and 7 fields an effective Hamiltonian con- 
sisting of bare A «-> ttN, ^yN vertices and energy-independent meson-exchange ttN potential 
(Fig. 3) and 'yN — > tcN transition operator (Fig. 5.). With the parameters listed in Table 1 
for the strong form factors and the bare mass of the A, the model can give a good description 
of irN scattering phase shifts up to the A excitation energy region. The only adjustable 
parameters in the resulting pion photoproduction amplitude are the coupling strengths Ge 
of the electric E2 and Gm of the magnetic Ml transitions of the bare A <-> ^yN vertex and 
the less well-determined coupling constant g W NN of the u meson. We have shown that the 
best reproduction of the recent LEGS data of the photon-asymmetry ratios of the 7p — > 7T°p 
reaction depends sensitively on these three parameters and yield Gm = 1.9 ± 0.05 and 
Ge = 0.0 ± 0.025 within the range 7 < g^NN < 10.5. Within these ranges of parameters, 
the predicted differential cross sections and photon-asymmetry ratios are in an overall good 
agreement with the data of 7p — > n°p, 7J9 — > ir + p, and — > n~p reactions from 180 MeV to 
the A excitation region. The model however underestimates the 7iV — > ttN cross section at 
energies above the A region. This is expected since the constructed model does not include 
inelastic channels, such as 7rA, pN channels, which should start to play a significant role 
at energies above about 350 MeV. Including these channels could also be needed to resolve 
the difficulty in fitting P13 tcN phase shifts (Fig. 22). The constructed effective Hamiltonian 
is free of the nucleon renormlization problem and hence is suitable for nuclear many-body 
calculations. 

We have also analyzed the i^-matrix method which is commonly used to extract em- 
pirically the 7iV — > A transition amplitudes from the 7^ — > ttN data. It is found that 
the assumptions made in the .fT-matrix method are consistent with our meson-exchange 
dynamical model. Our average value of the E2/M1 ratio Rem = (-1.8 ± 0.9)% is close to 
(-1.07 ± 0.7)% of Ref. 0. The helicity amplitudes calculated from our bare 7iV — ■> A vertex 
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are in good agreement with the predictions of the constituent quark models (Table IV). The 
differences between these bare amplitudes and the empirical values extracted from the data 
by using the iT-matrix method are shown to be due to the non-resonant meson-exchange 
mechanisms. This suggests that the bare vertex interactions in our effective Hamiltonian 
can be identified with hadron models in which the irN and nirN "reaction" channels (both 7r 
and N are on their mass-shell) are excluded in the calculation of the N* excitation. Unfortu- 
nately we are not able to pin down the E2/M1 ratio of the bare 7 — > A vertex by considering 
the existing data of photon-asymmetry ratios and differential cross sections. More precise 
data of other spin observables are needed to make progress. This will be pursued when 
the data becomes available, along with the extension of our approach to investigate pion 
electroproduction. 

The unitary transformation method developed here can be extended to higher energy 
regions for investigating higher mass N* resonances. To proceed, we need to perform the 
unitary transformation up to second order in the coupling constants to account for the 2-7T 
production channels. The resulting scattering equations will be defined in a larger coupled 
channel space N* © ttN © 7A © itttN. This research program can be carried out in practice 
since the numerical methods for solving such a Faddeev-type coupled-channel equations 
(because of the presence of the three-body ttttN unitary cut) have been well developed [|53 
Our effort in this direction will be published elsewhere. 
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APPENDIX: DERIVATION OF irN POTENTIAL 



To see how the Feynman-amplitude-like expressions of Eqs. (3.12)-(3.16) are obtained 
in our approach, we give a detailed derivation of nN potential from the following familiar 
Lagrangian 

J'ttNN 



ip(x) lb Yr^(x) ■ d„(/)(x) (Al) 



m,, 

h.c. 



where ip(x), ip A (x), and <f>(x) are respectively the field operators for N, A, and n; T is 
the A — > N isospin transition operator, [h.c] means taking the hermitian conjugate of 
the first term. By applying the canonical quantization procedure (see section III about the 
problem concerning the A field), we can derive a Hamiltonian from Eq. (Al). The resulting 
Hamiltonian can be written as 

H = H + Hj (A2) 

with 

Ej = Hf + H? . (A3) 

where H P (H®) describes processes which can (cannot) take place in free space. Explicitly, 
we can write in second quantization form 

H -^i^-^m dpdp ' d ' k (A4) 



. fnNA I rn N I m A 



m n V E N (p) V E A (p>) p 



- u£, T a u f k^ At, b ii ai a 5(p + k - p') 



where a^ a , bt, and At are , respectively, the creation operators for ir, N, and A states, a is 
the pion isospin index, u and cj m are respectively the spinors of Dirac and Rarita-Schwinger 
fileds. Clearly, H p describes the A <-> ttN real processes (similar to Figs, la and lb with the 
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change a — > A) which can take place in free space. On the other hand, the virtual processes 
(similar to Figs, lc-lf) are due to the following Hamiltonian 

1 1 



dp dp dk 



(A5) 



x 



JivNN 


/ m 


1 rn j 


m w V 




En( P ')\ 



-|m p -/7 5 t q ftupbp bp[ - 5(p' - p- k)a %a + 5(p' -p+ k)a± c 
+ up, 75 r a fivpbpdt [ - S(p' + p - k) a %a + S(p' + p + k) a± Q 
+ %'75 fr a updp>bp -5(-p' -p-h)a^ a + 5(-p' -p + k)a± a 



+ V7s tyr a vpdp, dt - 8{-p' + p - fc) + + p + k) , 



— « 



ATA 



m 



-< up 1 , upk^ At, bpS{p' -p+ k)at a 



. E N (p>) V E A (p' 
+ u$ ■ vpk^ At,dt 5{p' +p + k) at ta 

■Zk»b+Ap6(p>-p-k)a^ a 



+ % 



■fr 



NA 



m 



m A 



Up' UJ r - 



m n y E N (p') y E A (p) 
+ vp> fjt kf, dp, Ap8(-p' -p - k)a^ a 

Note that the above equation includes an anti-nucleon spinor v, which is included to maintain 
the relativistic feature of the starting quantum field theory. 

To proceed, we need to derive the unitary transformation operator S. By the procedures 
outlined in section II, S is related to Hf. Hence, the actual task of deriving the nN potential 
is to evaluate Eq. (2.23) from the Hf and Hf defined above. Let us first focus on the first 
four 7riViV coupling terms of Eq. (A. 5). We need to consider 



\i> = btat |0 > 

1 P ka 1 



l/> = #4J°> 



P ka' 

The allowed intermediate states in Eq. (2.23) are \n >= btjO >, bp~ m a\, a , a^jO > for the first 
term, and In >= d + ^ bt, b% at at, ,10 >, d + ^ bt,bt 10 > for the other three terms involving 

' 1 —pn p p ka k a ' —Pm p p i 

the anti-nucleon component v. Substituting these intermediate states into Eq. (2.23) and 
performing straightforward operator algebra, we then obtain 
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(k'a,p'\Hj \ka,p) 



m N 



X Up' 



i=i 



u 



where 



= m N ^ )t'T a >up n up> n 7 b U*tJ- 

&N{Pn) * 



with p n — k + p — k' + p'\ 



p , ,lh^T a up m up m -f 5 ft'r a ,- 



E N {p) - 


E N (p n ) + E n (k) 




1 


E N (Pn) 


-E N {p')-E 7r (k')_ 




1 


E N {p) - 


- ^(pj - E n (k') 




1 



E N (p m ) - E N (p') + E N (k) 



with p m = p — k' — p' — k; 

M (3) = ™" 



1 



77i / JS^'V U-ft, ^-p n 75 # r a - 
&N\Pn) * 



(A6) 
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M< 4 > 



■-— — -7s }tr a v.p m v.p m ^ ft r a ,- 

&N{Pm) Z 



-E N (p') - E N (p m ) + E n (k) 
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E N {p) + E N (p m ) - E n (k>) 



by using the properties that 

m N 



1 



E N (p) p p 2E N (p) 
m N I 



E N (pf^ 2E N (p) 
one can easily show that for an arbitrary p 



m N + -f E N (p) -7-p 
- m N + 7 E N (p) -7-p 



(A7) 
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1 

" 2E N {p n ) 
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p -E N (p n ) po + E N (p 
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m N 
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in - m N 



PO - E N (p n ) 

where p n = (po,p n )- 

By using Eq. (A8), we can combine various propagators in Eq. (A6) to obtain 



i=i A 
+ 75 f r a - 



1 



1 



{1> '+ ¥) - m N (^+^')-'m N 
1 1 



+ 



75 fT a 

75 ft'r a > 



where p = (E N (p),p), k = (E^k),^). The above result looks remarkably simple. It re- 
sembles very much the usual Feynman amplitudes, except that the intermediate nucleon 
propagator is the average of two Dirac propagators for the momenta evaluated using the 
incoming or outgoing tcN momentum variables. 

The evaluation of the A terms is much more involved, but yields a similar form as given in 
Eqs. (3.15)-(3.16). Similar derivations can also be carried out to define the nN interactions 
, Eq. (3.14), due to the p meson. 
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TABLES 

TABLE I. The parameters of the irN models. The units are 1/Fermi for cutoff parameters A a 
and MeV for the bare A mass. Model-L and Model-H are obtained respectively from fitting ttN 
phase shifts up to 250 MeV and 400 MeV. 



Model 


f 2 

4-7T 


A-nNN 


SpNN9pTTTT 


Kp 


Ap 


fnNA 




m/s. 


Model-L 


0.08 


3.2551 


38.4329 


1.825 


6.2305 


2.049 


3.29 


1299.07 


Model-H 


0.08 


3.7447 


39.0499 


2.2176 


7.5569 


2.115 


3.381 


1318.52 



K<.nnrac«.<u-i-L VA .(A-)-( 3 ^p r V , F p (q) = L^S , F nNA (k) = (a#^V 



TABLE II. The calculated -zrA^ scattering lengths( in unit of Fermi) are compared with the 
values determined in Ref. [40] 



Model-L Model-H Koch-Pietarinen 

Sll 0.1588 0.1737 0.173 ± 0.003 

5 3 i -0.1191 -0.1198 -0.101 ± 0.004 

Pii -0.0976 -0.0864 -0.081 ± 0.002 

P 3 i -0.0509 -0.0478 -0.045 ± 0.002 

P 13 -0.0363 -0.0383 -0.030 ± 0.002 

P 33 0.2523 0.2797 0.214 ± 0.002 
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TABLE III. The parameters for the <yN -> 
Fig. 5. 


itN interactions defined in Eqs. 


(4.16)-(4.20) and 


fnNN, ^-ttNN 


- 


Table I 


fnNA, A-nNA 


- 


Table I 


9pNN = yj 9 pN N 9 pirn 


- 


Table I 


J^pNN = A p , K p 


- 


Table I 


9pTfy 


0.1027 


Ref.[45] 


9lott^/ 


0.3247 


Ref.[45] 


9uiNN 


7 - 10.5 


See text 


A-luNN = A p 


- 


Table I 


Kuj = 




See text 


G M (0) 


1.85 - 2.0 


See text 


G e (0) 


± 0.025 


See text 
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TABLE IV. The magnetic M x + and electric E x + amplitudes of A — > 7./V transition at W = 1236 
MeV. Rem = E x + /M 1 + . The amplitudes are in unit of 10~ 3 (GeV)" 1 / 2 . The numbers in the upper 
and lower rows for each case are respectively from using (g u NN -,Gm ,Ge) = (10.5, 1.85, + 0.025) 
and (7., 1.95, -0.025). 

M x + E x + Rem Average 



r A ^ 7 7V (Bare) 175 -2.28 -1.3% 



184 +2.28 +1.2% 



r A -> 7 iv (Dressed) 257 -6.9 -2.7% 



258 -2.26 -O.c 



(0.0 ± 1.3)% 



(-1.8 ± 0.9)% 



TABLE V. Helicity amplitudes of the A — > 7^ transition at W = 1236 MeV are compared 
with the values from Particle Data Group (PDG) [45] and the predictions of constituent quark 
models of Refs. [47,48]. The amplitudes are in unit of 10~ 3 (GeV)^ 1 / 2 . The numbers in the upper 
and lower rows for each case are respectively from using (g^NN ,Gm,Ge) = (10.5, 1.85, + 0.025) 
and (7., 1.95, -0.025). 

A PDG Dressed Bare Ref. [48] Ref. [47] 

A 3/2 -257 ± 8 -228 -153 -157 -186 

-225 -158 

A l/2 -141 ± 5 -118 -84 -91 -108 

-126 -96 
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FIGURES 

FIG. 1. Graphical representation of the interaction Hamiltonians Hf of Eq. (2.7) and Hf of 
Eq. (2.8). 

FIG. 2. Graphical representation of the effective interaction Hamiltonians Hj P of Eq. (2.19) 
and Hj Q of Eq. (2.20). 

FIG. 3. Graphical representation of the interactions of the effective Hamiltonian Eq. (3.8) in 
the coupled 7riV © A space. 

FIG. 4. Graphical representation of scattering equations defined by Eqs. (3.23)-(3.30). 

FIG. 5. Graphical representation of the effective interactions Ta^n and v^m of Eq. (4.13). 

FIG. 6. Graphical representation of pion production amplitudes defined by Eqs. (4.22)-(4.24). 

FIG. 7. The irN phase shifts calculated from the irN Model-L of Table I are compared with 
the empirical values of the analyses of Ref . (3{| (open squares) and Ref . |38j (solid squares) . Tl is 
the pion laboratory energy. 

FIG. 8. The on-shell matrix elements of irN potentials defined by Eqs. (3.11)-(3.16). Tl is the 
pion laboratory energy. The notations are Np : vn d , Ne ■ vn e , p ■ v p , Ad : v& D , Ae ■ va e , and 
TOT is the sum. 

FIG. 9. The upper half shows the mass of the A state defined by Eq. (3.30). m& = 1299.07 
MeV is the bare mass of Model-L of Table I, Mr = 1236 is the experimental resonance position, 
and W is the irN invariant mass. The lower half shows the bare and dressed A — * ttN form factors 
defined by Eqs. (5.1) and (5.2) with f w NA = 1-3/^jva- k is the pion momentum in the A rest 
frame. 
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FIG. 10. The ttN phase shifts (5) and scattering amplitude (t(ko, ko, W)) with the irN invariant 
mass W = E^(ko) + E 7r {ko) in the P\\ channel. The solid curves are from the full calculation. The 
dashed curves are from the calculation including only the nucleon pole potential vn d of Fig. (3c). 

FIG. 11. The photon-asymmetry ratios i? 7 = da\\/da_L of jp — ► 7r°p at 90° degrees. E 1 
is the photon energy in the laboratory frame. The results in the upper half are from using 
Gm = 1-85, Ge = 0.025 and three values of g^NN- The results in the lower half are from us- 
ing Gm = 1-85, QujNN = 10.5 and four values of Ge- The data are from Ref. M]42], 



FIG. 12. The region of the parameters (g^NN ,Gm,Ge) for describing the data of 7iV — > irN 
reactions. See text for the explanations. 

FIG. 13. The photon-asymmetry ratios(a) and differential cross sections(b) for ^yp — > ir°p reac- 
tion at four angles in the center of mass frame. The solid and dotted curves are respectively fom 
the calculation using the parameters (g^NN, Gm, Ge) = (10.5, 1.85, +0.025), and (7., 1.95, —0.025). 



The data are from Refs. 1 41 . 42 ] . is the photon energy in the laboratory frame. 



FIG. 14. Same as Fig. 13, except for the jp — » ir + n reaction. The data are from Ref. [42] 



FIG. 15. The photon-asymmetry ratios(a) and differential cross sections(b) for the jp — > tt°p 
reaction at four angles. The solid and dotted curves are respectively from the calculations using 
the parameters (g^NN, Gm, Ge) = (10.5, 2.0, +0.025), and (7., 2.10, —0.025). The data are from 



Ref. [ 41 , 42 1 . E^ is the photon energy in the laboratory frame. 



FIG. 16. Same as Fig. 15, except for the jp — ► ir + n reaction. The data are from Ref. [42]. 



FIG. 17. Differential cross sections of 7p — ► ir°p, (a) *yp — > 7r + n, (b) and ^yn — > ir~p (c) 
reactions. The solid (dotted) curves are calculated by using Gm = 1.85(2.0). Both calculations 
using the same (g^NN, Ge) = (10.5, +0.025). The data are from Ref. Q|. E 1 is the photon energy 
in the laboratory frame. 
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FIG. 18. The predicted multipole amplitudes M x + and E±+ in the total isospin / = 3/2 channel 
are compared with the empirical values of Refs. |H| (solid squares) and ]57| (open squares). The 
parameters used in this calculation are Gm = 1-85, g^NN = 10.5, with Ge = 0.025 (solid curve) 
and -0.025 (dotted curve). E~ is the photon energy in the laboratory frame. 

FIG. 19. The M\+ and multipole amplitudes of the dressed vertex T 7 at^a defined by Eq. 
(4.24) and the bare vertex r 7 jv^A are compared. The dressed vertex is a complex function written 
as r(a) = |r(a)|e l< ^ a ) with a = M 1 +,£' 1 +. W is the jN invariant mass. 

FIG. 20. The predicted K-matrix defined by Eq. (4.28). The solid curves are the full calcula- 
tions. The dotted curves are from the resonant term. The dashed curves (denoted as B) are the 
contributions from the nonresonant term /c 77r defined by Eq. (4.29). W is the jN invariant mass. 

FIG. 21. The M]+ residues A of the if-matrix (Eq. (4.32)) calculated from the dressed A — > 
defined by Eq. (4.30) and the bare vertex are compared in the upper half. Their corresponding 

E 4. 

ratios Rem = jj^- are compared in the lower half. W is the 7^ invariant mass. 

FIG. 22. The 7riV phase shifts calculated from Model-L (solid curves) and Model-H (dotted 
curves) are compared. The data are from Refs. |38| (solid squares) and |?9) (open squares). Tl is 
pion energy in the laboratory frame. 

FIG. 23. Photon-asymmetry ratios and the differential cross sections for the jp — > tt°p (a and 
b) and jp — ► 7r + n (c and d) reactions. The solid (dotted) curves are from calculations using 7riV 
Model-L (Model-H). The parameters are (g^NN, Gm, Ge) = (10.5, 1.85, +0.025). The data are 
from Ref. 0Q. E 1 is the photon energy in laboratory frame. 
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